Six-electron semiconductor double quantum dot qubits 
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We consider a double-quantum-dot (DQD) qubit which contains six electrons instead of the usual 
one or two. In this spin qubit, quantum information is encoded in a low-lying singlet-triplet space 
much as in the case of a two-electron DQD qubit. We find that initialization, manipulation, and read- 
out can be performed similarly to the two-electron case, and that energy gaps remain large enough 
that these operations can be performed robustly. We consider DQD potentials with parameters 
chosen to be representative of current experimental capabilities. Results are obtained using two 
complementary full configuration interaction methods. 



The gate model for quantum computation requires the 
coupling of many individual quantum bits in order to per- 
form quantum "gates" which in turn implement useful 
quantum algorithmSfi"— There have been many proposed 
ways of realizing a single quantum bit (qubit), in a variety 
of material systems. Examples of experimentally demon- 
strated solid state qubits include the spin of a donor- 
bound electron^, the charge state of a single electron con- 
fined to a double-quantum-dot (DQD}^i^, and the total 
spin of two electrons confined to a DQDiii^ In this work, 
we consider an extension of this latter type: a qubit de- 
fined by the spin of six electrons confined in a DQD. This 
system has recently been realized experimentallyi^ 

A qubit can be defined by the low energy singlet and 
unpolarized triplet states of two electrons confined in a 
double-quantum-dot (DQD) potentialii The splitting be- 
tween singlet and triplet states, called the exchange en- 
ergy, determines the rate of qubit rotation about the z- 
axis of the Bloch sphere and can be controlled quickly via 
electric gates by "tilting" the DQD, that is, raising the 
potential of one dot relative to the other <^ Qubit rotation 
about the x-axis of the Bloch sphere occurs when there is 
an effective magnetic field gradient across the two dots, 
either naturally occurring (as in GaAs) or applied by 
micro-magnetsiS or g-factor engineeringii (as would be 
needed in isotopically pure silicon). DQD qubits can be 
coupled electrostatically to form entangling gates equiv- 
alent to the standard CPHASE or CNOT gates up to 
single qubit rotations^ Coherent single qubit rotations 
have been demonstrated in gallium arsenide (GaAs)^ 
as well as silicon-germanium (SiGe)^^ in lateral DQDs, 
and an entangling gate between qubits has been demon- 
strated in GaAsJ^ 

While the solid state environment has great promise 
for future scalability, semiconductor heterostructures are 
inherently "dirty" , containing many defects and impu- 
rities even using state-of-the-art fabrication techniques. 
Particularly problematic are charge impurities, which are 
present in concentrations of order 10"'^'' cm~^ in GaAs and 
10^^ cm~^ in Silicon MOS.— These impurities result in 
a potential landscape that is very rough when there are 
few conduction band electrons present, making it difficult 



to empty a quantum dot down to its last electron. Ad- 
ditionally, if the impurities are not static, the resulting 
fluctuating charge can cause qubit decoherenceJ^ While 
singlet-triplet qubits are more robust to charge noise than 
charge qubits, to the extent that the singlet and triplet 
states defining the qubit space have differing charge dis- 
tributions (e.g., during an exchange gate), charge fluctu- 
ations will still result in gate errors. 

A natural way to reduce the impact of randomly placed 
charge impurities is to consider quantum dots containing 
more than one electron. It has been shown that the ex- 
change energy of double quantum dots containing more 
than one electron per dot is less sensitive to a charge im- 
purity than that of a two-electron DQDjii In this work, 
we demonstrate that a six-electron DQD forms a viable 
qubit which can be initialized, manipulated, and read out 
in a similar manner to a two-electron DQD qubit. We 
show that in addition to the reduced sensitivity to charge 
impurities such many-electron qubits offer richer control 
capabilities than two-electron qubits. This furthermore 
obviates the need to deplete a dot of all but a single elec- 
tron, a process that at best requires substantial time and 
resources, and at worst is exhausting to the point that the 
device is deemed unusable. The choice of six electrons in 
a DQD is motivated by the fact that three electrons in 
a single dot result in a single valence electron (after fill- 
ing the lowest S-type shell) so that behavior similar to a 
DQD with one electron per dot is expected. (This argu- 
ment assumes a single valley, an assumption we operate 
under throughout this work.) Picturing a quantum dot 
as an artificial atom containing shells of electrons relies 
on the dot's orbital energy spacing being large compared 
to electron-electron correlation energy. This regime is ac- 
cessible to experiment, as evidenced by predictable shell- 
filling data in both GaA a^^d^ and silicon dots^, but is 
not necessarily the regime of all dots. In particular, the 
dot confinement along with the electron effective mass set 
the orbital energy scale, which need not be larger than 
the energy scale of the Coiflomb interaction, especially in 
the case of large dots. 

We investigate the spectrum of an idealized DQD po- 
tential using two complementary configuration interac- 



tion (CI) techniques. The CI method was first applied to 
multi-electron quantum dots within the quantum com- 
puting context over a decade ago,— with the qualitative 
conclusion that dots with an odd number of electrons 
should suffice as qubits. We first consider an unreal- 
istically small DQD for which both of our CI methods 
readily converge and which clearly illustrates the com- 
mon qualitative behavior found for all the cases of inter- 
est. Qualities of the six-electron spectrum related to the 
initialization, manipulation, and read-out are compared 
with those of a two-electron DQD which has the same 
valence electron tunneling rate, which we define below. 
The use of this metric helps to isolate the multi-electron 
effects from those which occur by simply changing the 
number of electrons in the system. We later calculate 
the six-electron spectra using parameters more compat- 
ible with today's experimentally realizable DQDs, and 
find that comparable qualitative behavior can be found 
with these larger-sized dots. We use GaAs material pa- 
rameters throughout the text. 
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FIG. 1: DQD potential along x-axis, V^(2;,0), given by Eq. [T] 
with L — 30 nm and Eo = 3.0 meV. V is the minimum of 
parabolas centered at x = ±L with curvature proportional to 
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I. SYSTEM 

The double quantum dot systems we consider are 
treated within the effective mass approximation with an 
electrostatic potential V given as the minimum of two 
parabolic potentials. 
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There is assumed to be strong vertical (z-direction) con- 
finement due to the semiconductor heterostructure, and 
we model only the lateral {x and y) dimensions. The po- 
tential of Eq.[l]is parametrized by e, the energy difference 
or "tilt" between the left and right minima, L, half the 
distance between the minima, and Eq ~ hu, the confine- 
ment energy of the dots (in the lateral directions) . The 
dots we consider are isotropic, that is, the saine confine- 
ment exists in both lateral directions. We discuss this 
restriction below in the context of applying our results 
to realistic dots. We have in mind e as the primary con- 
trol parameter, which is adjusted by varying the voltages 
applied to electrostatic gates. Changing e changes the 
exchange energy and thus the rate of the qubit rotation. 
Plots of the potential along the a;-axis for different e are 
shown in Fig. [TJ 

Typical dots in GaAs have Eq = 1 — lOmcV, L = 
25 — 100 nm. and e in the milli-electronVolt range. 



II. METHODS 

The configuration interaction algorithni^^ is well suited 
to modeling the multi-electron physics of small systems 
and has proven to be the method of choice when studying 



few electron systems in which multi-particle coherences 
are import an t ^^'^^i^'^" —. It consists of building antisym- 
metric multi-electron states from products of single elec- 
tron states (Slater determinants) and diagonalizing the 
Hamiltonian in this multi-electron basis. Since the di- 
mension of the multi-electron state space grows rapidly 
with the number of electrons TV and the size of the sin- 
gle particle basis M (roughly as mCn), even for a six- 
electron system the multi-electron basis must be "trun- 
cated" , that is, reduced from it's full dimension. The 
degree of accuracy achieved by a CI implementation is 
dependent on the single particle basis chosen (its size 
and form) and upon the truncation strategy and amount. 
Configuration interaction techniques, even when trun- 
cated, are variational^ meaning that they will never un- 
derestimate the ground state energy (or the lowest en- 
ergy state of a given symmetry class, e.g. triplet states). 
This property provides a means of comparing different 
CI implementations (see Appendix [B]) in the absence of 
a known exact solution, as the CI which produces the 
lower energy provides a stricter (better) bound on the 
ground state energy. 

In this work, the CI algorithm is used to solve a n- 
electron effective mass Hamiltonian of the form 
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where k is an effective dielectric constant, and the single 
particle Hamiltonian for the z*^ particle is given by 
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We have used f and p to denote the position and mo- 
mentum, respectively, of the i**^ electron and V is the 
single-particle potential (Eq. [1]) . The effective mass m* 
is generally a tensor but is assumed to be a scalar here. 



as it is in GaAs and even silicon when dealing with only 
the lateral directions. The magnetic field B = W x A 
is determined by the vector potential A. Wc take A in 
the symmetric gauge, A = B/2{—yx + xy). Coupling 
between spin and magnetic field is given by the effective 
g- factor g* multiplied by the Bohr magneton /xb. GaAs 
material parameters, m* = 0.067 me, g* = —0.44, and 
K = 12.9 are taken. 

We use two complementary CI implementations to 
study the six- and two-electron DQDs of this paper. The 
first, which we call CTl, uses the Fock-Darwin states as 
a single particle basis. The Fock-Darwin states are the 
exact single-electron eigenstates of a 2D parabolic dot in 
the presence of a constant magnetic field, and have the 
form 
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where Iq = Ib/{1/^ + ^7^2)1/4^ ^^ ^ ^hc/eB, uj^ = 
eB/m*c, L'^{x) is an associated Laguerre polynomial, 
and the associated single particle energies are 
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The second CI implementation, which we denote CI-2, 
uses a set of s-type Gaussian orbitals as a single particle 
basis. We denote the number of spatial Gaussian func- 
tions riG, which is half the number of elements in the 
single particle basis due to spin. Each spatial part of a 
Gaussian orbital is centered at some point (a;o,2/o) ^nd 
made gauge invariant by placing a magnetic-field depen- 
dent factor in its exponent, with the resulting form 



g{x,y) = 7V'e-"(^-^«) Q-o'iy-voY 



(6) 



where A/" is a normalization factor and B is the magnetic 
field. The vectors ro = (xo,yo) and a set the position 
and size of the Gaussian function, respectively. The sizes 
of the Gaussians are set to match that of the ground 
state of a single dot with the confinement energy Eq, so 
a = {m*Eo)/2T? . The positions are set on a regular rect- 
angular grid with spacing equal to one half of the Gaus- 
sian "radius" l/-\/a (a somewhat arbitrary value found 
to give good results in practice). 

Apart from the basis being used, the two methods are 
identical (though they were developed separately and are 
completely independent). Both truncate the full space 
of antisymmetrized product functions by keeping only 
those whose non-interacting energy (found by evaluating 
only the first term of Eq. [2]) is at most a specified cut- 
off energy, Ec, above the non- interacting energy of the 



ground state. We test each implementation for conver- 
gence by checking that the resulting multi-electron en- 
ergies are converged with respect to E^. This conver- 
gence is specific to the set of DQD parameters, and is 
faster for smaller dots since greater confinement energy 
means that the low-energy many-electron states are al- 
most entirely spanned by products of the lower orbital 
states kept by the CI. Figure [5] shows the convergence of 
the lowest single level for each method for the parame- 
ters Eq = 24.26 meV, L = 20nm, e = 0, and B = 2.8 T, 
which we motivate in section Hill below . Both CI methods 
converge to nearly the same energy, with the convergence 
of CTl showing an abrupt drop at Ec = 20fvjJc/'2.. This 
drop is due to the inclusion of the first excited orbital 
(Fock-Darwin state) with the same magnetic quantum 
number as each dot's valence orbital. The relatively large 
Coulomb matrix elements between these orbitals results 
in the observed reduction in ground state energy. The 
convergence of CT2 is more gradual since Gaussians are 
not the excited eigenstates of a dot, but we see a similar 
decrease in the ground state energy near Ec = 2{)fujJc/2 
which indicates the Gaussian basis is doing a decent job 
of spanning the excited single-dot levels. We find that 
the drop in energy near 20fiwc/2 is larger in the case of 
CI-1, and that it takes a much larger energy cutoff and a 
large number of Gaussians (26 in Fig. [5]) to obtain as low 
a ground state energy using CT2. If CI-2 is used with 
fewer Gaussians a drop near 20Swc/2 is still seen, but is 
smaller than for cases with a larger Gaussian basis. 

Since CI-1 uses the exact eigenstates of a single dot 
as its basis, it performs best when the confinement en- 
ergy of dots are large (small dots spatially) so that there 
is a large gap between the valence orbital and those be- 
neath it. This results in multi-electron wave functions 
which have a character similar to the valence electron or- 
bitals, as in the case of atomic physics. The second CI 
implementation (CI-2) performs comparatively less well 
at tight confinement, measured by the number of config- 
urations that need to be kept to achieve similar ground 
state energies. The Gaussian basis of CT2, however, is 
better able to produce multi-electron wave functions that 
arc much larger (spatially) than the single-electron solu- 
tions of the same potential. We use CI-1 to model cases 
where the confinement energy (i?o) is large, and CI-2 in 
cases where Eq is small. What constitutes a "large" or 
"small" Eq is determined by comparing the resulting en- 
ergies and convergence properties of both CI implementa- 
tions for that particular case. Both CI implementations 
were necessary for the results we present in this paper, 
and for each set of parameters below we indicate which 
implementation was used. 



III. 6-ELECTRON DQD QUBIT SPECTRA 

The configuration interaction method is applied to sev- 
eral DQD and single-dot systems to assess the viability 
of a 6-electron DQD qubit. In this section, we first ex- 



306 " 

^ 304 
> 

E, 302 
ta 
300 



ground 
first excited 



> 



bq 



308 
306 
304 
302 
300 
298 



5 10 15 20 25 
Ec iiicJc/2) 



ground 
first excited 



30 



10 15 20 25 30 35 40 45 50 
E^ (h(oJ2) 



FIG. 2: Convergence of the lowest two ^^ = levels as a 
function of the cutoff energy using CI-1 (upper) and CI-2 
(lower). In the lower plot, CI-2 uses a basis of jig ~ 26 
Gaussian functions, arranged on each dot in a 3 by 3 array 
with one additional element along each compass direction. 



plain and motivate the parameters that are used and then 
present the results. 



A. Parameter Choices 

A double quantum dot in our model is parameterized 
by a confinement energy Eq, an inter-dot separation 2L, 
the magnetic field i?, and a detuning or tilt e. We have 
in mind a DQD qubit where e is tuned by varying the 
voltages of electrostatic gates to control qubit rotation 
while the other parameters arc held fixed, and thus the 
majority of the figures below plot energies vs. e. We con- 
sider three different sized 6-e DQDs, given by Eq = 24.26, 
7.28, and 1.82 meV. These values oi Eq correspond to dot 
diameters of 13, 25, and 50 nm, respectively, where the 
diameter d is taken from the size of the harmonic oscilla- 
tor ground state d = 2yjh? /{mEo). While the diameter 
is a good indicator of the size of a single electron dot, the 
multi-electron dots wc study here may have larger spa- 
tial extent, and it is meant only to give a rough idea of 
the size of the confining potential. Experimental DQDs 
more typically confine electrons into puddles of around 
25 — 50 nm across, and so we understand the larger two 
dot sizes to be more realistic than the smallest size. In 
many cases we lock the vertical magnetic field B to -Eq 
by setting uic = 2aj/\/99. This relationship is chosen out 



of convenience for the numerical methods, to simultane- 
ously avoid near-degeneracy effects (which would occur 
at too small a magnetic field) and the overlap of differ- 
ent shells in energy (which would occur at large magnetic 
field) ^ 

We choose the inter-dot separation to keep the dots 
separated enough so that charge sectors remain distinct 
yet close enough so that there is appreciable tunneling 
between the dots. Specifically, we define the "valence 
tunneling" of a DQD with N electrons to be the energy 
difference between the ground and first excited state of 
the same system with A^ — 1 electrons. This estimates 
the energy required to move an electron from the valence 
level of one dot to the valence level of the other, and 
thus gives a metric of the tunneling between the dots. 
The valence tunneling falls off with increasing inter-dot 
separation at large L, and we choose L for each of the Eq 
values above by somewhat arbitrarily choosing a value 
of the valence tunneling which lies on the tail of the va- 
lence tunneling vs. L curve. For Eq = 24.26, 7.28, and 
1.82 meV wc choose values 0.2, 0.1, and 0.1 for the va- 
lence tunneling, which result in values of L equal to 20, 
35 and 60 nm, respectively (see Fig.|3l). (The decrease of 
valence tunneling as L is reduced smaller than the dot 
diameter is spurious and signals the breakdown of the CI 
method.) 

When wc compare 6-electron DQD results with 2- 
electron DQD results, our 2-electron calculations use the 
same values of Eq and B but vary L so that the va- 
lence tunneling remains the same as in the 6-electron 
calculation. Essentially, the DQD potentials used for 
our 2-electron calculations have more closely spaced dots 
to maintain the same wave function overlap as in the 
6-electron case (since the 2-electron wave functions are 
smaller than the 6-electron wave functions). The reason 
we fix the valence tunneling instead of using the same 
potential for 2 and 6 electrons is to mimic the situation 
in a real device: we suppose that when a realistic DQD 
is tuned between 2 and 6 electrons that the inter-dot 
tunneling is also tuned so as to obtain similar behavior 
of the DQD (e.g., charge stability lines and honeycomb 
charge sector structure). This results in values of L ap- 
proximately equal to 16.5, 28.5, and 47.5 nm for the three 
values of Eq above (see Fig[31). 

In summary, the DQD parameters just discussed are 
given in Table HI which we refer to as our "small-dot" , 
"medium-dot" , and "large-dot" DQD parameters as in- 
dicated. We use these sets of parameters, and variations 
around them, for all of the results to follow. 



DQD Parameters 



Name 



Small-dot 

Medium-dot 

Large-dot 



Eo 



24.26 meV 
7.28 meV 
1.82 meV 



L {Ge-) L {2e-) 



20 nm 
35 nm 
60 nm 



16.5 nm 
28.5 nm 
47.5 nm 



B 



2.8T 
847 mT 
212 mT 



TABLE I: Summary of double quantum dot parameters. 




FIG. 3: Valence tunneling as a function of half the inter-dot 
separation L, for Eq = 24.26, 7.28, and 1.82 meV (shown in 
(a), (b), and (c), respectively). A vertical magnetic field B is 
set such such that hwc = 2JfJo/\/99 {uJc ~ eB/m). Horizon- 
tal lines indicate the constant value of the valence tunneling 
chosen to define the L parameters in each case. Calculations 
were performed using CI-2 with no ~ 18 and Ec — 55 meV. 



Only un-polarized, i.e., Sz = states are shown, and we 
arbitrarily cut off the upper part of the spectrum. We 
say a state is in "charge sector" (n, m) when there are 
n electrons in the left dot and m electrons in the right 
dot. The four slopes seen in Fig. |4] correspond to the 
four possible charge sectors: (3,3), (2,4), (1,5), and (0,6). 
(There are always equal or more electrons in the right dot 
since positive tilt (e) corresponds to the right dot having 
lower energy.) We denote the 5"^ = energy levels by 
a capital S, T, or Q, indicating a spin-0 singlet, a spin- 
1 triplet state, or a spin-2 quintuplet state, followed by 
a zero-based index indicating the energy ordering. For 
example, the lowest lying singlet and triplet states are 
denoted SO and TO, respectively. The progression of the 
lowest singlet eigenstate (SO) among the charge sectors is 
noted in Fig. SI The relative "size" of the regions where 
SO is in each charge sector is an indication of the shell 
structure. For example, the tilt required to remove an 
electron from a 3-electron dot and add it to another 3- 
electron dot is less then the tilt required to remove an 
electron from a 2-clcctron dot and add it to a 4-clcctron 
dot. 




FIG. 4: 6 electron energies vs. detuning using CI-1 for small- 
dot parameters (Eq = 24.26 meV, L = 20nm, B = 2.8 T) and 
energy cutoff Ec — ?>Qf>Wc/2 = 73.15 meV. 



In the small-dot DQD the orbital and multi-electron 
state energy spacings arc large enough to result in a clear 
and relatively simple spectrum which we present first. 
Then, since the small-dot parameters are unrealistic for 
actual devices, we focus on the medium- and large-dot 
DQD systems and consider how the many-electron spec- 
trum evolves as one moves into a parameter regime which 
more closely reflects real devices. 



B. Small-Dot DQD Spectrum 

The energy spectrum of a 6-electron DQD using our 
small-dot parameters (see Table |l]) is shown in Fig. S) 



We find that for our small-dot parameters the lowest 
two 5^ = states are a singlet and triplet state, and that 
there is a clear order-meV gap to higher states. When 
the singlet state of one charge sector crosses the singlet 
state of another there is an anti-crossing with a finite gap 
as shown by Fig. [51 which zooms in on each of the singlet 
anti-crossings. Each anti-crossing is qualitatively similar 
to the (l,l)-(0,2) crossing in the corresponding 2-electron 
spectrum shown in Fig. [5] Thus, the most basic criteria 
for a qubit are satisfied: there exists an isolated singlet- 
triplet ground space with variable splitting throughout a 
wide range of tuning. We will discuss qubit character- 
ization further below, but first let us consider how the 
energy spectrum changes with respect to the dot size. 




FIG. 5: 6 electron energies vs detuning using CI-1 for small- 
dot parameters and energy cutoff Ec = 30huJc/2 = 73.15 meV, 
zoomed to each of the three anti-crossings, (3,3)-(2,4), (2,4)- 
(1,5), and (l,5)-(0,6) of the lowest singlet state, "SO". 




non-zero magnetic field along z. Thus, for qubit oper- 
ation analogous to the 2-electron DQD case, we require 
a singlet ground state in addition to the existence of an 
isolated 2-level manifold with a tunable energy splitting. 




15 20 25 
e (meV) 

FIG. 7: 6 electron energies vs detuning using CI-2 {no — 18, 
Ec = 55 meV) for medium-dot parameters {Eo = 7.28 meV, 
L = 50 nm, B = 847 mT). ^z = for aU states, and labels 
denote spin (S=singlet, T=triplet) and energy-ordering index. 



FIG. 6: 2 electron energies vs detuning using CI-1 for small- 
dot parameters {Eo = 24.26 meV, L = 16.5 nm, B = 2.8 T) 
and energy cutoff Ec = 30fta;c/2 — 73.15 meV. The left panel 
shows many of the 2-electron energies over a large range of 
detuning; the right plot shows a zoomed view of the lowest 
singlet anti-crossing. S(l,l) and S(0,2) label the lowest singlet 
state in the (1,1) and (0,2) charge sectors, respectively, and 
(T_,To,T+) label the three lowest (1,1) triplet states with 
Sz = —1, 0, and 1, respectively. 



C. Medium- and Large-dot DQD Spectra 

Having seen that the 6-electron DQD can be operated 
as a qubit in the usual way in the case of small dots, we 
now turn to the question of whether this is also true of 
DQDs with larger dots. This amounts to a question of 
energy scales, since by enlarging the dots we are decreas- 
ing the energy spacing between the orbitals of a single dot 
and consequently the spacing between the many-electron 
levels. Ideally the spectrum would have an isolated 2- 
level Sz = ground space formed by a singlet and triplet, 
or by a singlet and the unpolarized state of a higher-spin 
multiplet, with the singlet state at lower energy. The rea- 
son wc desire a singlet ground state is for qubit initializa- 
tion: in order to initialize a Sz = state there must be a 
regime accessible to the qubit in which a singlet state is 
the true ground state with respect to all Sz sectors. If a 
higher-spin state is the ground state, e.g. a triplet state, 
then initialization is made more challenging if not impos- 
sible because it will always be energetically preferable to 
initialize into one of the polarized triplet states for any 
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FIG. 8: 6 electron energies vs detuning using CI-2 {no ~ 18, 
Ec — 55meV) for medium-dot parameters, zoomed to each 
of the three anti-crossings, (3,3)-(2,4), (2,4)-(l,5), and (1,5)- 
(0,6) of the lowest singlet state. Sz = for all states, and 
labels denote spin (S=singlet, T=triplet) and energy-ordering 
index. 



Previously we considered dots with the confinement 
energy Eq = 24.26 mcV, which corresponds to a ground 
state electron wave function (in GaAs) with diameter ap- 
proximately 13 nm. In the medium-dot and large-dot pa- 
rameter cases wc now consider, the ground state wave 
function diameters are approximately 25 nm and 50 nm 
respectively. The energy spectrum for the medium-dot 
DQD parameters is shown in Fig. [71 with zoomed-in 
views of the singlet anti-crossings in Fig. |8l The spec- 
trum is very similar to the small-dot case. There exists 
a low-lying isolated singlet-triplet subspace with avoided 
crossings, with an mcV-scale gap to excited singlet and 




FIG. 9: 6 electron energies vs detuning using CI-2 {no ~ 10, 
no cutoff) for large-dot parameters (Eq = 1.82 meV, L — 
50 nm, B — 212 mT). 3^ = for all states, and labels denote 
spin (S=singlet, T=triplet) and energy-ordering index. 
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FIG. 10: 6 electron energies vs detuning using CI-2 {nc ~ 10, 
no cutoff) for large-dot parameters, zoomed to each of the 
three anti-crossings, (3,3)-(2,4), (2,4)-(l,5), and (l,5)-(0,6) of 
the lowest singlet state. Sz = for all states, and labels 
denote spin (S=singlet, T=triplet) and energy-ordering index. 



triplet states. The singlet-triplet gap between ground 
states in the (2,4) region is small and negative, but can 
be tuned to a large positive value by the magnetic field, 
as will be discussed in Section IV below. The (0,6) ex- 
change energy gap is positive and sufficiently large to 
enable proper qubit initialization. 

Moving to even larger dots with Eq = 1.82 meV brings 
us into a regime where the CI technique is no longer 
fully convergent. In the previous two cases where Eq = 
24.26 meV and 7.28 meV, the CI is able to converge with 
respect to its truncation, and in the Eq = 7.28 mcV case 
we get qualitative convergence with respect to the Gaus- 
sian basis size (see Appendix |A|). In particular, adding 
more Gaussian elements may lower the state energies a 
small amount but does not cause level reordering. In 
the present case of large-dot DQDs (Eq ~ 1.82 meV), 
this is not the case, and special care must be taken in 
the interpretation of the results. The numerical chal- 



lenge stems from the fact that a large number of single- 
particle states are required to capture even the lowest 
energy multi-particle states, so that less truncation can 
be applied and we are therefore limited to smaller Gaus- 
sian basis sizes. This limitation, in turn, restricts the 
variational freedom of the many-particle states, leading 
to less certainty in qualitative features such as the rela- 
tive spacing and ordering of energy levels. In the present 
case, computational power limits us to a Gaussian basis 
of size uq = 10 if we require complete convergence with 
respect to truncation. 
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FIG. 11: Energy spectrum of a single dot with 6 electrons. 
Dot confinement energy Eq — 1.82 meV. At a vertical mag- 
netic field B — 212 mT, the spectrum approximates the (0,6) 
regime of a DQD with our large-dot parameters. Computed 
using CI-2 with no ~ 13 and no energy cutoff. Sz = for 
all states, and labels denote spin (S=singlet, T=triplet) and 
energy-ordering index. 

The resulting spectrum is shown in Fig. [HI with mag- 
nifications near the anticrossings shown in Fig. [TOl As 
expected, we see the continued trend of smaller energy 
spacings in the low-energy spectrum. We also note that 
a non-singlet state is the lowest-lying state in the (2,4), 
(1,5), and (0,6) charge sectors for nc = 10, and that such 
a situation presents a potential problem for initialization 
as discussed above. A non-singlct ground state is plausi- 
ble since decreased orbital energy separations can make 
it favorable for electrons to occupy higher orbitals and 
form non-singlets in order to reduce Coulomb repulsion. 
This line of reasoning is verified deep in the (0,6) region 
of the spectrum, where we can approximate the DQD 
system as a single dot with six electrons and solve just a 
single dot problem using the CI. This allows us to place 
more Gaussians in the single dot and ultimately gives us 
a converged CI spectrum, which is shown in Fig. 1111 We 
see that at i? = 212 mT a triplet state is the ground state. 
However, it is important to note that if the magnetic field 
were reduced below 200 mT or increased above 500 mT, 
the ground state will become a singlet, so that the abil- 
ity to initialize the qubit can be restored by tuning the 
magnetic field appropriately. 

Another observation that can be made regarding the 



DQD spectrum shown in Fig. [10] is that there may be 
intervening levels between the two qubit states. For in- 
stance, in the (2,4) region the lowest three states are all 
triplets, while the lowest singlet state is the third ex- 
cited state. Thus if we were to use the lowest singlet and 
triplet states as the qubit space, there would be two in- 
tervening (triplet) levels in the (2,4) region. This would 
make it difficult to operate the qubit in the (2,4) region 
(a possibility the details and benefits of which will be 
elaborated on in Section IV) , unless it is possible to tune 
the magnetic field such that the qubit states form a well 
separated ground space in (2,4) while maintaining the 
singlet as the ground state in (0,6). If this is not possi- 
ble, then one could still restrict qubit operation to the 
(3,3) region, tuning through (2,4) solely for the purposes 
of initialization and readout. 



IV. QUBIT CHARACTERIZATION AND 
CONTROL 



dot DQD parameters is shown in Fig. [121 The qualita- 
tive shape of this curve is generic for 6e- DQDs with 
well separated dots: there are two plateaus correspond- 
ing to the regions where both singlet and triplet are in 
the (2,4) or (0,6) charge sector, interspersed by two flat 
"valleys" in the (3,3) and (1,5) regions. The heights of 
the (2,4) and (0,6) plateaus are related to the single-dot 
exchange energies for dots with 2, 4, and/or 6 electrons, 
and the sharpness of the ascent to and descent from the 
plateaus are a function of the inter-dot tunneling ampli- 
tude. Larger tunneling (a lower barrier) results in greater 
spacing between the singlet and triplet avoided crossings 
at a given charge sector transition, and thus a more grad- 
ual ascent/descent. The relatively small exchange ener- 
gies of the (3,3) and (1,5) regions are determined by the 
inter-dot coupling. These features are analogous to the 
2e case apart from the fact that there are two plateaus 
instead of one (as there are four charge sectors instead 
of two) . The possible advantages due to the existence of 
two plateaus will be discussed below. 



A. Similarities between 6- and 2-electron DQDs 

The control of a 6-electron singlet-triplet qubit is 
nearly identical to the 2-electron case. By tilting the 
DQD, the exchange energy J (the difference between the 
lowest triplet and singlet energies) can be tuned to cause 
qubit z-rotation, and an existing or applied gradient in 
the magnetic field causes an x-rotation. Initialization is 
performed by biasing the qubit so that there is a large 
gap between the singlet and triplet states and waiting for 
the DQD to relax to the ground state. Read-out is done 
by biasing the qubit so that the singlet and triplet are 
in different charge sectors and measuring the charge sec- 
tor via a nearby charge sensor (usually a quantum point 
contact or a single-electron-transistor). 
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FIG. 13: Exchange energy of a single dot with confinement 
energy Eq in magnetic field B = 847 mT (the same value used 
in our set of medium-dot DQD parameters, see Table[l]) for 2, 
4, and 6 electrons. Computed using CI-2 with no ~ 13 and 
no energy cutoff. 
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FIG. 12: Comparison of the exchange energy vs detuning in 
6- and 2-electron DQDs using small-dot parameters (given in 
TableU]). Note that the dot-separation parameter L is different 
for the two cases in order to equate the valence tunneling of 
the 6- and 2-electron DQDs as described in Sec. lIIIR . We find 
comparable orders of magnitude, indicating that z-rotations 
of this 6-electron DQD qubit can be controlled in a similar 
manner to the rotation in the corresponding 2-electron DQD 
qubit. Calculations were performed with CI-1 using Ec = 
30&.i;c/2 = 73.15meV 



The exchange energy as a function of tilt for our small- 



We also compare with a two-electron result in Fig. [T^l 
The two-electron DQD has a different DQD potential 
(L = 16.5 nm instead of 20 nm). in order to match 
the valence tunneling of the six-electron DQD as de- 
scribed above. We see that the magnitude of the six- 
electron DQD exchange energies is similar to that of the 
two-electron DQD. Although it is computationally time- 
consuming to repeat the calculations across the whole 
range of dot sizes, it is easy to show that the similar- 
ity between six-electron and two-electron DQD exchange 
plateaus is independent of dot size by plotting the ex- 
change energy of a single dot as a function of its confine- 
ment energy. When the inter-dot coupling is small, the 
heights of the double-dot exchange plateaus are essen- 
tially the exchange energies of a single-dot at the appro- 
priate filling, so that the 2-electron (0,2) and 6-clcctron 
(0,6) plateau heights correspond to the two- and six- 
electron exchange energies of a single dot. Figure [T3l thus 
implies the similarity of 6-electron (0,6) and 2-clcctron 
(0,2) DQD exchange plateaus over a large range of con- 
finement energies. We can also estimate the dependence 



of the (2,4) plateau height based on dot confinement en- 
ergy when (1) the 2- and 4-electron single dots each have 
an isolated singlet-triplet low energy subspace, and have 
substantially different exchange energies, and (2) the dots 
are weakly coupled. In this particular case, the lowest 
6-electron singlet or triplet states can be approximated 
by a singlet in the dot with the larger exchange energy 
and a singlet or triplet, respectively, in the dot with the 
smaller exchange energy. For B = 847 niT. the case of 
Fig. [131 the 2-electron exchange energy is only somewhat 
larger than the 4-electron exchange energy (which is neg- 
ative due to the degeneracy of the first excited dot or- 
bital ~ Hund's rule), and we see that while the (2,4) 6- 
electron exchange energy is also negative, it's magnitude 
is much smaller than the 4-electron single-dot exchange 
energy, indicating the 6-electron singlet and triplet have 
excited state components beyond the simple picture out- 
lined above. 
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FIG. 14: Comparison of the singlet-triplet energy coupling, 
h, in 6-electron (top pane) and 2-electron (bottom pane) 
DQDs using small-dot parameters (given in Table [J). We find 
comparable order of magnitudes, indicating that x-rotations 
of this 6-electron DQD qubit can be controlled in a simi- 
lar manner to the rotation in the corresponding 2-electron 
DQD qubit. Calculations were performed with CI-1 using 
Ec = 30hujc/2 = 73.15 meV. 



In order to design qubit operations, it is very help- 
ful to construct an effective Hainiltonian describing the 
two-level ground space. In the case of our multi-electron 
DQDs, this Hamiltonian assumes the following form: 



H, 



eff 



J 

2 



h 



(7) 



where a^^ o^ are Pauli matrices, J is the exchange 
splitting between the lowest two levels \g) and |e), and 
a magnetic field gradient with Zeeman energy hjl = 
{g\ ^'''^r^^ x\e) couples these states, with AB = Br-Bl 
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FIG. 15: Comparison of the singlet-triplet energy coupling 
h, and exchange energy J, vs detuning in the (3,3) region 
of a 6-electron DQD (top pane) and the (1,1) region of the 
corresponding 2-electron DQD (bottom pane) using small-dot 
parameters (given in Table [l]). The similarity indicates that 
this 6-electron DQD could be controlled as a qubit in the same 
way as the corresponding 2-electron DQD. 



denoting the difference between the magnetic fields in the 
right and left dots. In the case of two-electron singlet- 
triplet qubits, field gradients on the order of tens to hun- 
dreds of milliTesla have been generated and stabilized 
via a feedback mechanism involving repeated driving 
through singlet-polarized triplet anticrossing a^^'^^ . Us- 
ing CI-1, we have computed the Zeeman energy as a 
function of detuning for both the two-electron and six- 
electron cases with AS — 60 mT, as shown in Fig. [HI 
It is apparent from Fig. [15] that the six-electron values of 
J and h in the (3,3) region are comparable to those in 
the (1,1) region of the two-electron case (both are on the 
order of ^eV). Therefore, the effective two- level Hamil- 
tonians for two and six electrons are qualitatively the 
same, so that qubit rotations can be implemented using 
similar procedures in the two cases. For example, recent 
progress in developing robust singlet-triplet qubit rota- 
tions in two-electron DQD o^^'^^ should apply equally well 
to the six-electron case. 

Qubit read-out is performed via a nearby charge sen- 
sor, and relies on the sensitivity of the charge sensor to 
the DQD state. Many factors enter into the coupling 
between the DQD charge state and the signal from a 
charge sensor, and we do not account for many of these 
here. Instead we consider the difference in Coulomb en- 
ergy of the singlet and triplet states at a point in space 
50 nm away from a DQD (this distance is arbitrary, and 
is meant as an estimate of the location of a charge sens- 
ing constriction or dot). We define Coulomb energy to 
be the integral /e^|^(r)p/(r — ro)dr, where ^ is the 
wave function whose Coulomb energy we're computing 
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FIG. 16: Comparison of the difference in Coulomb energy 
(defined as /e'^|^(r)p/(r — ro)dr) between the singlet and 
triplet states 50nm away fi"om the DQD for the 6-electron 
(top) and 2-electron (bottom). Small-dot parameters are used 
as given in Table|I]and computation was done with CI-1 using 
Ec = 30&Jc/2 = 73.15 meV. 



(the lowest singlet or triplet state) and r^ is the point in 
space 50 nm away from the DQD. Thus, the difference in 
the Coulomb energy of the singlet and triplet states at 
the location of a charge sensor gives an indication of how 
strongly the DQD state couples to the sensor (or equiv- 
alently, how easy it is to perform read-out on the qubit). 
We find that the ease of read-out is similar for 2-electron 
and 6-electron DQDs (using our small-dot parameters), 
shown in Fig. 1161 

Clearly, in order to use the (2,4) or (0,6) regions for 
readout or initialization requires adiabatically sweeping 
through avoided crossings in the spectrum (see, e.g.. 
Fig. [8]). By the Landau- Zener formula for an avoided 
crossing with minimum energy gap 25, the rate at which 
the energy gap is changed, ■;;, should be less than 6'^ /h 
to avoid the diabatic transition that occurs with prob- 
ability e~'^'^^ ^^^ . Thus, an adiabatic sweep across an 
anticrossing starting from a detuning with a larger en- 
ergy gap 2 A requires ramping time T > /iA/(5^. Even 
for A '-^ 2(5 ^ 100/xeV, ramp times of tens of ps satisfy 
this condition. On the other hand, if the ramp is too slow 
the magnetic field gradient will rotate the singlet toward 
the triplet during the sweep. For typical field gradients 
less than a ^eV, this prohibits sweeps longer than about 
a ns. So, sweeping the detuning for readout and initial- 
ization should not pose a problem for ramping times on 
the order of 100 ps. In the event it does become difficult 
to satisfy the above conditions, nonlinear ramps provide 
an obvious alternative approach. 



B. Differences between 6- and 2-electron DQDs 

Along with the many similarities just discussed, there 
are some important differences between 6-electron and 
2-electron singlet-triplet DQD qubits. One, of course, is 
the difference that motivates this work: more electrons 
result in an increased ability to screen charge impurities, 
as investigated in Ref. llTl The incremental increase in 
screening ability with the number of filled shells is found 
to fall off at high electron number, but there is clear 
benefit from the increase to 6-electrons (1 filled shell) 
from 2-electrons (no filled shells) in a DQD. This may 
mean that 6-electron DQD qubits are less susceptible to 
charge noise and the random potential caused by disor- 
der, making them easier to realize from both a materials 
and charge-noise standpoint. Additionally, if the num- 
ber of electrons could be inferred without emptying the 
DQD completely and counting up charge transitions (as 
it is typically done at present), a 6-electron DQD could 
prove easier to fabricate from a lithography standpoint. 

On the other hand, although the energy scales are ac- 
ceptable in the 6-electron DQDs, the low-lying qubit sub- 
space is less isolated from the higher levels than in the 
2-electron case, particularly when dealing with large dots. 
This is the only disadvantage we see to using 6-electron 
DQD qubits, and we believe it is far outweighed by the 
benefits afforded. Clearly, though, one cannot increase ad 
infinitum, the number of filled shells to achieve increased 
robustness to charge impurities. The many-body energy 
gaps are dependent on the dot size (confinement energy) 
and magnetic field in addition to the electron number, 
and what energy gaps can be tolerated will depend on 
temperature, qubit gate speed, and readout technique. 
In short, precisely determining the optimal (odd) num- 
ber of electrons per dot will be a non-universal system- 
dependent quantity. With that said, we speculate that 
for typical GaAs and Silicon (assuming a large valley 
splitting) DQDs this optimal number will be somewhere 
between three and nine electrons per dot, since three ap- 
pears to be clearly better than one, and by seven the 
charge-screening benefits are already seen to be minor 
(c.f. fourteen-electron DQD results using the frozen-core 
approximation in Ref. llTI ) . 

Another difference that is potentially very useful is 
that a 6-electron qubit has more charge sectors and hence 
more control capabilities than a 2-electron qubit. There 
is a practical question of how much tuning range can 
be realized in an actual 6-electron device, but for the 
present discussion let us assume that tuning between 
all four charge sectors is possible. Then there are two 
plateaus where J is relatively small (typical parame- 
ters give J ranging from nano- to micro-electron-volts), 
and two plateaus with potentially much larger J (milli- 
electron- volts with typical parameters) . As we discussed 
above, for DQDs with well separated dots, the heights of 
these plateaus are of the same order as the exchange ener- 
gies of single dots with appropriate numbers of electrons. 
The plateau heights depend on the confinement energy. 
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as shown in Fig. [131 and the magnetic field. The i?-ficld 
dependencies of the (2,4) and (0,6) plateaus are not pro- 
portional to each other, so the ratio of the plateau values 
can be tuned by B-field, as shown in Fig. [T7]at medium- 
dot DQD parameters. The existence of regions where 
the exchange energy is largely insensitive to the detun- 
ing of the DQD can be used to rotate (or not rotate) the 
qubit in a way that is robust to detuning noise. This is 
the long-standing idea behind the "sweet spot" or "park" 
state, and has been studied^^i^"— for the case of rotation 
in a 2-electron DQD. It was found that due to current 
limitations of controlling electronics (minimum voltage 
pulse widths), the value of the exchange flat must fall in 
the range of 1 — 10 ^eV, which is much smaller than the 
0{ meV) exchange energy required for good initialization 
and read-outj^ In the 6-elcctron case, one can tunc the 
system so that the one plateau is low (1 — lO/icV) and 
one is high 0{ meV) in magnitude so that both operations 
can be done at a single magnetic field (see, e.g., Fig-ITT]). 
Another new capability with 6 electrons is to use both 
(2,4) and (0,6) plateaus for rotations at different rates, 
or even different signs (again, sec Fig. [IT)) . This could be 
leveraged to perform error correction to the extent that 
the environments seen by the (2,4) and (0,6) states arc 
similar—. 
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FIG. 17: Exchange energy vs detuning using CI-2 {nc = 10) 
for medium-dot parameters Eq = 7.28 nieV and L = 50 nm, 
and magnetic fields B = 500, 847, and 1000 mT. The ability 
to tune the two exchange "plateaus" in the (2,4) and (0,6) 
regions with magnetic field is one possible advantage of a 6- 
electron DQD qubit over its 2-electron counterpart. 



V. CONCLUSIONS 

We have shown that filling a DQD with 6-clcctrons 
results in a qubit that can be just as easy to initial- 
ize, manipulate, and read out as the well-established 2- 
electron qubit. Large exchange energies in the (2,4) and 
(0,6) charge sector regions allow for initialization, and the 
exchange energy and singlet-triplet coupling via i?-field 



gradient [h] dependencies arc similar in the two cases, 
indicating qubit manipulation proceeds in the same way. 
We find similar read-out sensitivity, since in both cases 
within the "readout regions" located between charge sec- 
tors there is a single electron charge transfer that differ- 
entiates the singlet and triplet states. We find exchange 
energies and gaps to excited non-qubit-spacc states are of 
order tenths of milli-electron Volts ( 3> fcT in a dilution 
refrigerator) in DQDs sized to match experimental capa- 
bilities. Furthermore, 6-electron DQDs offer an enhanced 
charge screening as well as richer control capability. This 
work implies that the devotion of experimental resources 
to the realization of single-electron dots is unnecessary 
for, and even detrimental to, qubit performance. In 
summary, comparison of 6-electron and 2-electron DQDs 
indicates that multi-electron quantum dots may be su- 
perior than single-electron dots for quantum computing 
purposes due to their additional controlability and intrin- 
sic screening of unwanted charge impurities. 
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Appendix A: Convergence of CI-1 and CI-2 

We consider in this appendix the convergence of the 
CI-1 method as a function of the cutoff energy, and of 
the CI-2 method as a function of the number of Gaus- 
sians used and the cutoff energy. First, we consider the 
convergence behavior of 6-electron DQDs when using the 
small-dot parameters in Table HI The convergence of CT 
1 with cutoff was shown above (Fig. ^, and consists of 
steps at the energies where higher-shell Fock-Darwin or- 
bital are introduced into the basis. At Ec = 30huJc/2, the 
value used to produce Figs. S] and [5l the single particle 
basis included the lowest ten states of each dot, consist- 
ing of Fock-Darwin S, P, D, and F states, and the many- 
body basis consisted of 1636 states. In the case of CI-2, 
Fig. [TH] shows the convergence of the ground state en- 
ergy within the (3,3) and (0,6) regions of detuning space. 
The flattening of the curves at high Ec indicates conver- 
gence with respect to the cutoff used. Note, however, 
that even after results at a given basis size (fixed nc) 
are converged with respect to the cutoff, increasing the 
number of Gaussian basis elements no can continue to 
lower the energy of the ground state. The disparity be- 
tween the _Bc-converged values at different no is greater 
in the (0,6) detuning regime than in the (3,3). This is 
a confirmation of what one may intuitively expect: that 
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when all six electrons are in a single dot the additional 
variational degrees of freedom given by increasing the 
number of Gaussians per dot has more energetic impact. 
In the (2,4) and (1,5) regimes we see convergence behav- 
ior intermediate between the (3,3) and (0,6) regions, and 
thus focus our convergence analysis to the (3,3) and (0,6) 
regimes. 
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FIG. 18: Convergence of the lowest singlet energy (SO) of the 
CI-2 method for a 6-electron DQD in the (1,1) regime (upper 
pane) and (0,6) regime (lower pane). Small dot parameters 
Eo = 24.26 meV, L = 20 nm, and B = 2.8 T, are used and the 
tilt e is indicated above each plot. 
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FIG. 19: Lowest singlet level vs detuning of a 6-electron DQD 
using CI-2 with different numbers of Gaussian basis elements 
no- Small dot parameters are used: Eo — 24.26 meV, L — 
20 nm, B — 2.8 T, and for each value of no the energy is 
converged with respect to the cutoff parameter E^. In the 
10-Gaussian case, the Gaussian's centers are arranged 5 per 
dot in a "plus sign" pattern, with one Gaussian in the center 
of the dot and one along each of the -(-/-x and +/-y axis 
directions. In the 18 Gaussian case, the centers are arranged 
in 3 by 3 grids on each dot, and in the 26-Gaussian case in 
addition to a 3 by 3 grid there is a Gaussian even further out 
along each of the -f /-x and +/-y directions. 



Next, we turn to the convergence behavior of the 
"medium" and "large" -dot DQDs of the text, given by 
the 6-electron DQD parameters of Table H] Figure HOI 
shows the convergence of the CI-1 ground state at zero 
tilt. It is difficult to conclude that convergence is ob- 
tained from these plots, as possibility of small steps be- 
yond Ec ~ SOhuc/^ cannot be eliminated, and we prefer 
the CT2 for these dot sizes because of its smoother con- 
vergence and lower ground state energy. 



It initially appears concerning that the ground state 
energy is not converged with respect to the size of the 
CI's single particle basis, that is, the number of Gaus- 
sians. While ideally one would like to sec convergence in 
uq as well as Ec, we fine that this impossible from a com- 
putational standpoint and that the lack of convergence 
does not affect the qualitative aspects of the results. Fig- 
urc ll9l providcs an example of this, showing the 6-clectron 
spectrum for uq = 10, 18, and 26 each converged with re- 
spect to Ec- We see that all three curves show the same 
qualitative behavior over the entire range of detuning, 
with the larger no curves lying slightly lower in energy. 
Thus, while we require results to be converged with re- 
spect to the CI's cutoff energy we do not (and cannot, 
given practical computational limitations) require that 
the results are numerically converged with respect to the 
single particle basis size. Instead, we content ourselves in 
a qualitative convergence with respect to the basis size 
and the layout of the Gaussians. We remark that in the 
CI-1 method, the basis size is also fixed to a pre-defined 
number of Fock-Darwin orbitals, and that we make the 
same decision to prioritize convergence with respect to 
Ec in that case, even though this is not as explicit in the 
text. 
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FIG. 20: Convergence of the lowest two energies of the CI- 
1 method for a 6-electron DQD using the "medium" (upper 
pane) and "large" (lower pane) dot parameters of TableUwith 
zero tilt (e = 0). 



Figures [3T] and [H] show the convergence behavior of 
CI-2 in the (3,3) and (0,6) regions of detuning similar to 
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Fig. [THl We notice progressively less step-like jumps in 
the energy curves as the dots become larger, indicating 
that there is the shell-structure picture is becoming less 
valid and relevant single particle levels are more numer- 
ous and spread out over energy. We again note the larger 
differences in converged energy in the (0,6) regime, and 
that in the large-dot DQD case the energy if ng = 18 is 
barely just barely converging with respect to Ec- This 
motivates the care taken in the text in the interpretation 
of the large dot results. 
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basis elements. Figure [23l compares the energy of the low- 
est singlet and triplet states (SO and TO) when CI-2 uses 
10 Gaussians and 26 Gaussians. The agreement is worst 
in the large tilt region, and CI-1 finds a lower ground 
state energy by a few milli-clcctron volts in this regime 
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FIG. 22: Convergence of the lowest singlet energy (SO) of the 
CI-2 method for a 6-electron DQD in the (1,1) regime (upper 
pane) and (0,6) regime (lower pane). Large dot parameters 
Eo — 1.82 meV, L — 60 nm, and B = 212 mT, are used and 
the tilt e is indicated above each plot. 



FIG. 21: Convergence of the lowest singlet energy (SO) of the 
CI-2 method for a 6-electron DQD in the (1,1) regime (upper 
pane) and (0,6) regime (lower pane). Medium dot parameters 
Eo = 7.28 meV, L = 35 nm, and B = 847 mT, are used and 
the tilt e is indicated above each plot. 



Appendix B: Comparison of Fock-Darw^in and s-type 
Gaussian CI bases 

In this appendix we compare the two configuration in- 
teraction implementations, CI-1 and CI-2, used in this 
work. While not of direct relevance to the results other 
than giving further confirmation and corroboration of 
each implementations correctness, we feel a detailed com- 
parison will still be of interest in its own right. Compar- 
isons are made using the "small-dot" DQD parameters of 
the text, since for these parameters both CI-1, which uses 
Fock-Darwin states as a single-particle basis, and CI-2, 
which uses a array of s-type Gaussian functions as a basis, 
converge well with cutoff energy Ec and the number of 



in the 10-Gaussian case. In the 26-Gaussian case the 
agreement is good to a small fraction of a milli-electron 
volt throughout the detuning range. The 26-Gaussian 
CI-2, however, is much more expensive than the CI-1 in 
terms of the number of configurations considered and to- 
tal computation time. This highlights the fact that the 
CI-1 is better suited to small-dot DQDs as discussed in 
the text. 

For larger dot sizes, such as the medium and large dot 
parameters in the text, we find that the CI-2 method 
produces lower ground state energies than CI-1. We be- 
lieve this is due to the fact that products of Fock-Darwin 
orbitals are less good variational approximations of the 
true multi-electron wave functions in larger dots, whereas 
the single particle basis used by CI-2 is more generically 
appropriate. 

Overall, the agreement between the two independently 
developed CI codes gives us confidence in the correctness 
of the implementations and lends credence to the accu- 
racy of a converged result using either implementation. 
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FIG. 23: Comparison of 6 electron energies produced by the 
two different CI implementations, CI-1 and CI-2, when CI-2 
uses no = 10 s-type Gaussians in the upper plot and no = 26 
in the lower plot. Small-dot parameters [Eq = 24.26 meV, 
L = 20nm, B = 2.8 T) are used. 
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